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Abstract 

The timing and dynamics of many diverse behaviors of mammals, e.g., patterns of animal foraging or human 
communication in social networks exhibit complex self-similar properties reproducible over multiple time scales. In this 
paper, we analyze spontaneous locomotor activity of healthy individuals recorded in two different conditions: during a 
week of regular sleep and a week of chronic partial sleep deprivation. After separating activity from rest with a pre-defined 
activity threshold, we have detected distinct statistical features of duration times of these two states. The cumulative 
distributions of activity periods follow a stretched exponential shape, and remain similar for both control and sleep deprived 
individuals. In contrast, rest periods, which follow power-law statistics over two orders of magnitude, have significantly 
distinct distributions for these two groups and the difference emerges already after the first night of shortened sleep. We 
have found steeper distributions for sleep deprived individuals, which indicates fewer long rest periods and more turbulent 
behavior. This separation of power-law exponents is the main result of our investigations, and might constitute an objective 
measure demonstrating the severity of sleep deprivation and the effects of sleep disorders. 
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Introduction 

Sleep deprivation 

Although good sleep (like nutrition and physical exercise) is 
considered a basic contributor to human health and well-being, its 
chronic deprivation seems symptomatic of modern societies. In the 
last century, the average sleep duration has shortened from 
approximately nine hours to no more than seven hours in many 
countries, e.g., USA [1], UK [2], or Japan [3]. Some recent 
research shows that both total sleep deprivation and chronic sleep 
reduction may lead to similar effects in terms of physiological, 
affective and cognitive consequences (e.g., [4,5]). The most 
obvious effect of sleep loss is the daytime drowsiness, an 
underestimated problem concerning operators of transportation 
and in other "critical-safety" work settings. Apart from drowsiness, 
sleep deprivation involves impaired immune, endocrine and 
metabolic functions, and profound neurocognitive deficits. A list 
of cognitive sleep loss consequences set in a review by Durmer and 



Dinges [6] includes a variety of symptoms: from slowed reactions, 
omission and commission errors, and a decline in working 
memory performance, to deterioration in divergent thinking and 
increased likelihood of unproductive problem-solving. Cognitive 
domains are, however, affected diversely by sleep loss, so that 
sustained attention deteriorates much more than the performance 
of challenging working memory tasks [4]. Impairments in 
performance are accompanied by changes in performance self- 
ratings [5]. 

Daytime drowsiness is characterized by an urge to sleep, a lack 
of energy, and decreased ability to complete tasks, often 
resembling the state of alcohol intoxication. Slowed movements, 
reduced facial expressions and muscle tone are typical physical 
symptoms of being sleepy. Thus, the spontaneous locomotor 
activity may be considered as an index of alertness and sleepiness. 
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Scale-free distributions 

Power laws are ubiquitous in nature and have been repeatedly 
detected not only in physics, biology, earth and planetary sciences, 
but also in economics and finance, demography, epidemiology, 
and social sciences [7-11]. They are typically observed in the 
vicinity of continuous phase transitions, where the underlying 
physical processes and fluctuations of measured physical observ- 
ables exhibit self-similarity at all scales. Scale invariance is 
therefore commonly considered a signature of "criticality" 
indicating the complexity of the system and slow decay of 
spontaneous fluctuations. Of special interest are the scale-free 
distributions in time, e.g., the waiting-time distributions of light 
and dark states in quantum dots [1 2] , or dwell-time distributions in 
closed biological ion channels [13], where the inverse power-laws 
of state-duration times have been observed and identified with 
Poisson shot noise in blinking (opening) events. Accumulating 
evidence also demonstrates [8,14—16] that the dynamics of 
spontaneous behavior exhibits scale invariance. Time recordings 
of locomotor activity in rodents and humans have shown that the 
spatial and temporal distribution pattern of fluctuations appear 
unchanged regardless of the time scale of observation, thus 
pointing to the aforementioned universal scaling laws. In 
particular, recent studies by Sun et al. [17] indicated that the 
scaling exponent of the power law detected in temporal 
autocorrelation of activity significantly correlates with the severity 
of Parkinson's disease symptoms. Similarly, universal scaling laws 
have been found in locomotor activity periods of humans sufiering 
from major depressive disorders [8]. The disruption of the 
characteristic universality classes of such laws has been further 
addressed by Proekt et al. [14] in studies on dynamics of rest and 
activity fluctuations in light and dark phases of the circadian cycle. 

This paper continues along this line of research, seeking to 
determine standards for measurable criteria by discriminating 
human behavioral organization as impaired by sleep deficiency. In 
particular, the aim of the present study is to find whether the same 
individuals show measurable differences when undertaking two 
difierent styles of everyday Kfe. The chosen experimental setup 
and subsequent statistical analysis allow us to observe such 
disparities in locomotor activity between subjects sleeping regu- 
larly and those undergoing chronic sleep deprivation. As discussed 
further, those distinct variations in patterns of activity can be 
already detected in the first day of sleep deficit, leading to 
significantly different scaling exponents for fluctuations of 
actigraphy recordings in both groups. When interpreting the 
findings we also point to their potential application in identifying 
patients with sleep disorders. 

Materials and Methods 

Actigraphy measurements were performed on healthy individ- 
uals over one week of their normal life [rested wakefulness (RW)] 
and one week of partial sleep deprivation (SD) (access to the data: 
[18]). The circadian cycle of both groups differs substantially: 
while RW individuals have relatively long "nights" and short 
"days", members of the SD group are characterized by a r(;vcrsed 
pattern of longer "days" and shorter "nights", which clearly 
influences their activity/rest patterns. To overcome this problem 
we normalized the "days" and "nights" of both groups to the same 
length, as explained in the Data analysis subsection. The 
resulting time series were statistically analyzed and compared to 
former studies [8] performed in the "normal" (RW) phase. Since 
bouts of activity/rest obey different distributions of duration, the 
best choice of a threshold(s) diSerentiating between the two states 



seems to be crucial, and is thoroughly discussed in the Results 
section. 

Participants 

Twenty four paid volunteers (12 females and 12 males; mean 
age 22.7 years, S.D. = 1.6) participated in the study. They were all 
healthy, non-smokers, and drug-free. They were asked to Hmit 
alcohol and caffeine intake during the experimental weeks. They 
reported regular sleep pattc'rns and no sleep-related problems, 
controlled with Pittsburgh Sleep Quality Index [19] and Epworth 
Sleepiness Scale [20]. Participants were informed about the 
procedure and goals of the study, and provided their written 
consent. The study was approved by the Bioethics Commission at 
Jagiellonian University. 

Among the 24 sets of data collected 1 7 were selected for further 
analysis (8 females and 9 males; mean age 22.8 years, S.D. — 1.8), 
the other 7 were corrupted by either removing a recording device 
during the experiment, or not following the sleep schedule. 

Data acquisition 

The data acquisition comprised of one week of unrestricted 
sleep according to individual needs, i.e., rested wakefulness (RW), 
and one week of daily partial sl(x:p d(:pri\ ation (SD), with a two- 
week gap in between the two measurements. Half of the subjects 
began with the RW phase followed by the SD phase, while the 
other hcdf had the order reversed. During the sleep deficit week, 
the participants were asked to shorten their sleep by 33% of their 
'ideal sleep' by delaying bed-time and using an alarm clock in the 
morning. The precise length of the restricted sleep was calculatc-d 
individually for each participant, where the indi\'idual sli;ep need 
was determined on the basis of the questionnaire administered 
before entering the experiment ("If you were totally free to plan 
your day and had no duties at all, at what time would you go to 
sleep and get up?"). In half of cases the self-reported length of sleep 
was verified with actigraphy before entering SD phase. 

For the 17 selected subjects the average unrestricted sleep, as 
measured by actigraphy, was 8 h 16 m + 44 m (mean + S.D.); in 
the sleep deprivation conditions it was curtailed by 2 h 20 m + 
56 m (27.7% ±9.6%) and amounted 5 h 57 m + 37 m. 

Movement tracking was recorded with Micro Motionlogger 
SleepWatch (Ambulatory Monitoring, Inc., Ardsley, NY), worn on 
the participant's non-dominant wrist. The data were collected in 
1 -minute epochs in the Zero-Crossing Method (ZCM) mode, 
which counts the number of times per epoch that the activity 
signal level crosses zero (or a threshold very close to zero). The 
working limitation of the ZCM mode is the difficulty in registering 
the acceleration of movements, which may potentially cause high 
frequency artifacts to be counted as a considerable movement. 

Data analysis 

Day and night or wakefulness and sleep. The raw 

actigraph data collected in the ZCM mode represent the number 
of counts per epoch as a function of time. The first step in the 
analysis is the localization of the 'day' ^ 'night' transition (or 
wakefulness <=? sleep). This is done by the procedure consisting of 
smoothing the raw data over longer periods of time (usually tens of 
minutes) and counting actigraph activity above a predefined 
threshold Tskep- The periods where rest is predominant (above 
given percentage hgieep) are counted as sleep. If they are 
interrupted by instantaneous activity, e.g., due to a change of 
body position, they are glued together. The resulting period is 
defined as sleep or 'night'. The remaining time is called 
wakefulness or 'day'. 
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Next, one has to avoid artifacts connected with forced, different 
lengths of sleep/wakefulness in the two week period under study. 
Whereas the RW individuals sleep on average 8 hours and are 
alert 1 6 hours per day, the SD ones sleep approx. 3 hours less and 
stay alert 3 hours longer. To make the two samples comparable, 
out of the total 24 hour period we take 5 hours of sleep, which are 
common to all participants, and 16 hours of consecutive 
wakefulness, again common to both groups. We call it the (5+ 
16) mode. For comparison, we also analyze the reversed 
combination: 16 hours of wakefulness followed by 5 hours of 
sleep, denoted as {16+5). 

Threshold(s) separating states of activity and rest. To 
precisely distinguish the rest periods from periods of activity, a 
threshold separating the two states must be chosen. We resolve this 
problem in two different ways: with a single or double threshold. 
In the single threshold method we select one specific value of ZCM 
activity, TzcM, which best separates the two types of behavior 
investigated. We also test the sensitivity of analyzed distributions to 
changes of TzcM around the optimal value. 

The double threshold procedure, in turn, is a hysteresis-like 
method, in which two step values Tjown and T^p are introduced. 
We then define the beginning of the activity period as the moment 
when ZCM activity exceeds Tup, and the end of this period when 
the activity falls below Tjown. At that moment the state of rest 
begins, and it terminates when surpassing T^p. We compare the 
results following from these two approaches. 

Statistical analysis. The raw actigraph data X{t), split into 
activity and rest, are subject to further statistical analysis. We 
count the number of activity/ rest periods of a given duration and 
calculate the resulting probability density function (PDF) p(t) of 
duration time T. To better assess the statistics of rare events in tails 
of PDF's we construct, as the main measure of the discussed 
phenomena, the (complementary) cumulative distribution C(«) of 
duration lengths a: 



C(a) = Pr(T> a) E 



p{x)dx. 



(1) 



The function represents the survival probability for the system 
to stay in a given state up to the time a. For a stationary time series 
the survival probability C{a) is expected to have a characteristic 
scale (relaxation time Xrel) related to the probability per unit time X 
to undergo a change of the state (see Appendix). Put differently, 
C(a) is then a simple exponential function of dwell times, 



C{a)-- 



with X=\ l%rel- 



In order to check the degree of temporal correlations in the 
activity recordings, the power spectrum of the signal was 
derived from the Fourier transform of the signal correlation 
function 



lim-< 



exp {ift:)X{}') exp (iff)X{t")dt'dt"}, (2) 



and evaluated for several time periods. The parameter / denotes 
the inverse of time f = t~^. 

The numerical estimates of cumulative distributions were fitted 
with two mathematical formulae: a power-law of the form 

Cia)~a-^ (3) 
for rest periods and a stretched exponential form 



C(a)~ exp( — aa^ 



(4) 



for activity periods. The fitting was performed using log-log or log- 
hnear data, respectively, in order to account for the tails in the 
distributions. The fitted parameters a, P, and y were then 
compared for several combinations of time periods and subjects. 

Each cumulative distribution was constructed from rest/ activity 
periods collected over a single day [i.e., a single (5+16)- or [16+5)- 
hour period] from all the 1 7 participants in RW or SD condition. 
The statistics for individual subjects is insufficient to construct 
reliable cumulative distributions. Consequently, for each day we 
obtained a single data point characterizing the distribution of the 
whole group in RW condition and a single point for SD condition. 

Throughout the paper the parameters characterizing a single 
cumulative distribution of rest/ activity periods are given together 
with the standard error of the fit (e.g., exponent t] of power 
spectra); whenever we take an average of the fitted parameters 
over the whole week, the weighted mean - marked with a 
horizontal bar - together with the standard deviation is given (e.g., 
y). In order to compare RW and SD groups, such means of 
parameters fitted for several consecutive days were compared with 
two-tailed Student's /-tests performed at 95% confidence level, 
preceded by a set of tests for equal variances. These tests were 
performed with Mathematica 9.0 software [21] (functions TTest 
and VarianceEquivalenceTest). 

Since the cumulative distributions, and consequentiy the 
parameters, may change depending on TzcM or T^p and Tjo„„, 
the thresholds were chosen so as to optimize goodness of fit of the 
two models (3)-(4) and thus best separate rest and activity periods. 
Whereas fitting the stretched exponential did not show any 
discernible optimum in the range TzcM = 50 — 150, fitting the 
power law had a clear optimum in several goodness of fit measures 
(namely, minima in sum of squares error Err, '•/} statistic, Akaike 
and Bayesian Information Criteria). As a result we chose the 
numerical values of thresholds 7zcM = 85 and 



= 85, Tdot 



= 60 based on this criterion. The details of the 



fitting procedures and measuring the goodness of fit are described 
in File SI. 

Results 

Raw data and their spectra 

The raw actigraph data XiJ) collected in the ZCM mode for an 
exemplary participant are shown in Fig. 1. The pattern of the 
activity events in RW mode (top-left panel) exhibits a highly-visible 
circadian rhythmicity, whereas episodes of activity measured in the 
SD mode (top-right panel) present a substantial portion of a short- 
interval chattering. The separation of "days" and "nights" has 
been performed with T^ieep = 40 and h^ieep = 0.8. The red solid and 
dashed lines above the data mark our normalized selection of 
periods under investigation, (5+16) hours and (16+5) hours, 
respectively. Additionally, in the bottom panels of Fig. 1 we plot 
the time series of increments Af = X{t+l) — X{t) normalized by 
its standard deviation (S.D.). To assess the nature of fluctuations in 
A(, we analyzed distributions of A//S.D. and the representative 
frequency histogram (which, by construction, stands for estimation 
of the probability density function of scaled fluctuations A,/S.D.) is 
drawn in the bottom panel of Fig. 2 along with a Gaussian 
probability density function of the same mean and variance. The 
"experimental PDF" is clearly leptokurtic: pronounced hea\'y tails 
show that by comparison to the Gaussian distribution, fewer low- 
amplitude fluctuations and many more large ones are observed in 



PLOS ONE I www.plosone.org 



3 



September 2014 | Volume 9 | Issue 9 | el 07542 



Scale-Free Fluctuations in Behavior of Sleep Deficient Humans 



Control period, subject no 24 

1_6+5h 

5+16h 



Sleep deprived, subject no 24 




2000 



4000 6000 

t [min.] 



8000 



10000 



4000 6000 

t [min.] 



10000 



Figure 1. Examples of typical activity recordings. (Left panel) Activity of a control subject (RW mode), and (right panel) of a sleep deprived 
subject (SD mode). The overall nonzero activity counts X{t) are depicted on the vertical axis. The lower row displays the time series of standardized 
increments Ai = X{r+l)-X(t). 
doi:1 0.1 371/journal.pone.01 07542.g001 



experimental sets, thus demonstrating their strong deviation from 
the standard central limit theorem. 

The degree of correlation between subsequent intensities of 
activity X(t) is well represented by its spectral density S(f). The 
illustrative plots for a single individual and different time span (one 
week and 24 hours) are drawn in Fig. 2. For the RW data sets, 
corresponding to a couple of hundred events, the power spectra 
exhibit a clear universal, algebraic scaling law S(f)ccf^'', with the 
scaling exponent rj close to 1 = 1.03 + 0.02 (left panel) and 
ri= 1.09 + 0.02 (right panel)]. The respective values for subject 24 
in the SD mode are f; = 0.99 + 0.02 (week) and (7 = 1.06 + 0.02 
(day). The middle row of the same figure displays the power 
spectra obtained for increments with slopes if xt] — 2 confirm- 
ing long-range correlations in the fluctuations of intensities X(t). 

Consequendy, these data demonstrate that the analyzed 
stochastic process X(t) representing variability in intensity of 
activities is not of the white-noise type, in which case the 
consecutive events would have been memory-free and the 
correlation function would have been a Dirac delta function with 
the corresponding spectral power S(f)ccf°. 

Cumulative distributions and their fits; single threshold 
analysis 

With the choice of the threshold Tzcm = ^^> the cumulative 
distributions C{a) of rest period durations, averaged over all 
participants, are plotted in Fig. 3a for RW and SD modes in the 
(5+16) combination. Analogous graphs in the (16+5) combination 
are depicted in Fig. 3c. Lidividual curves represent data for 
consecutive days. To minimize possible influence of individual 
variations among participants, we have replotted the same 
cumulative distributions C(a/a) as functions oi a/a, where a is 
the individual average, cf Fig. 3b and d. All curves stay close to 
straight lines on the log-log scale over 2 order of magnitude in 
time, which is a signature of power-law behavior of the form given 



by Eq.(3). We performed fits of the cumulative distributions for 
each curve and present the extracted exponent y in Fig. 4. There 
are several conclusions which can be inferred from this analysis. 
Firstly, we do not see significant differences between the (5+16) 
and (16+5) recordings. Secondly, the y values of the control group 
are close to those reported earlier [16] for healthy humans 
(y = 0.9 + 0.13 versus average y = 0.78 + 0.05 in this study). 

Finally - and we consider this observation our main result - the 
analysis performed indicates a significant difference in behavioral 
motifs between the control group and sleep deprived individuals. 
The higher coefficient y = 0.85 + 0.03 derived for sleep-deficient 
individuals emphasizes the fact that the pattern of their resting 
times consists of more short periods and, respectively, fewer longer 
inactivity time intervals than in the control group. This 
observation contrasts with the results for the rest-time distributions 
of depressed humans [8,16], where lower scaling exponent y, and 
thus heavier tails in the cumulative distribution, were observed for 
disordered individuals. The difference in the exponent y that we 
observe is significant: a two-tailed Student's f-test between sets of 
resulting y coefficients for the RW and SD groups was performed 
at 95% confidence level (p = 0.0015 and 0.030<0.05 for (5+16) 
and (16+5) settings, respectively) providing evidence for a 
statistically notable difference between means of these two groups. 
It is interesting that after the change of the scaling parameter y, 
already on the first day of sleep deprivation (Fig. 4), we do not 
observe any additional trends during the subsequent days. One 
should also notice that even if the exponents are close to each 
other on some days (e.g., days 5 and 7 in Fig. 4 right panel), the 
distributions are clearly distinguishable from each other, as can be 
seen in Fig. 5, where C{a) in the rest state for SD and RW 
individuals are plotted for exemplary days 6 and 7. 

In an analogous manner, the cumulative distributions of the 
activity periods for RW and SD cases are plotted as functions of a 
on log-log and log-lin scales in Fig. 6a,c; furthermore. Fig. 6b,d 
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A,/S.D. 



Figure 2. Temporal universality of signal X{t) depicted for a typical subject (No. 24). Top panels show the spectral densities S(f) 
evaluated (left panel) for the one-week experimental time series (i.e., top panels of Fig. 1 ) and (right panel) for 24 hr periods averaged over the week. 
In both cases slopes of fitted power laws (dashed lines) exhibit 1//'' behavior with the characteristic exponent f|= 1.03 + 0.02 (left panel) and 
f/ = 1 .09 + 0.02 (right panel), respectively. In the middle two panels the similar spectral analysis is shown for the time series of increments (i.e., bottom 
panels of Fig. 1) exhibiting 1//'' scaling with exponents rj' ^ —0.73 + 0.02 (left panel) and ij'= — 0.71 +0.02 (right panel). These results clearly indicate 
that the activity events are long-range correlated and the corresponding stochastic process of switching between active/non-active periods is not 
memoryless (IVIarkovian). The peaks located at 1/23.99 hours, marked by red triangles in the left panels correspond to the circadian rhythm. The 
bottom panel shows the frequency histogram of standardized increments A,/S.D. (black - unrestricted sleep; red - sleep deprivation) compared to a 
Gaussian probability density function of the same mean and variance (dashed line) on log-linear scale. 
doi:1 0.1 371 /journal.pone.01 07542.g002 



presents those distributions as functions of a/a. Both original and 
rescaled cumulative distributions of activity periods collapse well 
onto the stretched exponential form Eq.(4). The values of the fitted 



parameters a and P for RW and SD groups are not significantly 
different {p>0.05, two-tailed t-test). This result remains in 
agreement with previous studies on humans and rodents [8,15,16]. 
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Figure 3. Cumulative distribution of resting times in Kiuman motor activity. Double logarithmic plots represent (panels a, c) cumulative 
distribution C{a) of rest periods as a function of duration time a and (panels b, d) cumulative distribution C(a/a) of rescaled rest periods as a function 
of rescaled duration time a/a for all RW subjects (black symbols) and all SD (red symbols). Here, a stands for the individual mean of rest period 
duration. Each curve corresponds to one of six consecutive days. Continuous lines show slopes of the fitted power-laws. Top panels (54-16) mode; 
bottom panels (164-5)mode. The threshold separating state of activity versus rest rzcj/ = 85. 
doi:1 0.1 371/journal.pone.01 07542.g003 
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Figure 4. Exponent of the rest-periods distributions for RW (blaci<) and SD subjects (red) as a function of consecutive days. The 

cumulative distribution assumes (over almost two decades) a power-law form C{a)Ka^'''. Notation as in Fig. 3. In both cases the characteristic index y 

is significantly lower for SD. 

doi:1 0.1 371/journal.pone.01 07542.g004 
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Figure 5. Sample distributions of rest periods with fitted slopes. The continuous lines are best fits; their extent indicates the range of data 
points fitted (the first and last fifteen data points are not fitted). Panel (a) is an example where both the distributions and the fitted slopes can easily 
be distinguished. Panel (b) shows an example in which the fitted slopes are nearly equal, although the difference of distributions is clearly visible. 
doi:1 0.1 371 /journal.pone.01 07542.g005 



Activity periods, log-log scale 



Rescaled activity periods, log-log scale 



Control 

Sleep deprived 
5+16h 




1.000 
0.500 



0.100 
0.050 - 



0.010 
0.005 



0.001 - 



Control 

Sleep deprived 
5+16h 



(b) 




500 



Activity periods, log-linear scale 
— Control 



0.05 0.10 0.50 1.00 5.00 10.00 

a/a 

Rescaled activity periods, log-linear scale 



Sleep deprived 
5+1611 




T 1 .000 ^ 
0.500 



Control 

Sleep deprived 
5+16h 



100 



200 300 

a rmin.l 



400 



500 




a/a 



Figure 6. Cumulative distributions of activity periods. The curves on log-log and log-linear scales follow the typical pattern for a stretched 
exponential C(a)= exp( — W) as a function of (panels a, c) the length of activity periods a or (panels b, d) the length a rescaled by the individual 
average a. All curves collapse on a similar stretched exponential function; there is no significant difference between RW and SD subjects. The 
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The choice of thresholds separating the activity and rest periods 
is an important issue. By selecting one particular value of ZCM 
activity we separate two different temporal distribution profiles: 

the power law for states of rest and (stretched) exponential for 
states of activity. It is thus clear that there must be an optimal 
value that best distinguishes the two. The criterion we used to find 
it was the best fit of the resting-state cumulative distribution C(a) 
to the power-law form. The resulting sum of squares error Err, 
along with the goodness of fit expressed by x^, are plotted in Fig. 7 
as functions of the threshold value Tzcm- 

Robustness of the results: Cumulative distributions and 
their fits in double threshold analysis 

All the results discussed above were obtained with the single 
threshold definition (Fig. 1-7). In order to further test the 
sensitivity of the results to the definition of a threshold, we refined 
our analysis by considering a hysteresis-like distinction of the 
locomotor activity (see section Data analysis). With the choice 
Tdown = 60 and Tup = 85 we plot the cumulative distributions C{d) 
of rest (Fig. 8a,c) and activity (Fig. 8b,d) time intervals for both 
RW and SD groups, together with the fitted exponents of the 
power-law cumulative distribution of resting states (Fig. 9). In the 
case of activity periods no clear difference between the RW and 
SD individuals was observed. The overall fit (averaged over days 
and individuals) yields for the RW sample a = 0.31 + 0.04 and 
^ = 0.49 + 0.03. Similar values for the SD mode are not 
significantly different (p>0.05; two-tailed t-test). In contrast, the 
difference between the control and sleep deprived individuals can 
be well observed in the profiles of the rest-periods cumulative 
distributions. Numerically, the power-law fit yields y = 0.75 + 0.05 
for RW and y = 0.86 + 0.03 for SD samples. Altogether, we 
observe that cumulative distributions obtained with this definition 
agree with those obtained in the single thresholding method, 
therefore confirming the robustness of the results obtained. 

Discussion 

The evidence of motion recorded by actigraphs is broadly used 
in medical practice aimed at understanding typical circadian- 
rhythm patterns of healthy subjects, and at detecting possible 
rhythm disorders in humans. Although the method is not indicated 



for the routine diagnosis, the actigraphy recordings may be useful 
in quantitative evaluation of the variety of neuropsychiatric 
diseases. Specific analytical methods applied to actigraphy records 
can reflect various features of physiological activities and serve as 
an outcome measure in characterizing disturbances of the 
circadian-rhythm patterns and alterations in locomotor activity 
in specific populations like children (apnea testing) or individuals 
suffering from depression [16], Parkinson's disease [17], schizo- 
phrenia [22], or dementia [23]. 

In the former studies by Nakamura et al. [8, 16], actigraphy was 
employed to determine the difference in human behavioral 
organization between healthy subjects and patients with major 
depressive disorders. The authors analyzed activity and resting 
time durations and found that the cumulative distributions follow a 
stretched exponential form for activity periods and scale-free, 
power-law behavior for resting time intervals. Moreover, by 
analyzing the behavioral organization of humans suffering from 
major depressive disorders, they found significantiy lower power- 
law scaling exponents for the rest-period durations than for 
healthy control groups. 

In analogy to those investigations, we have performed the 
experiment and applied statistical evaluation of the results, keeping 
in mind that our RW and their healthy control groups should 
reach similar conclusions. Therefore, e.g., the choice of activity 
threshold which separates state of activity from rest has been 
validated by the use of a criterion similar to [8,16], i.e., the overall 
average of nonzero activity counts. 

Evaluation of residence-time distributions and calculation of a 
mean waiting time to an activity event are broadly studied 
problems in biological physics and have been addressed in a 
number of studies [9,13,14,24—29]. Stochastic systems with 
internal states in which discrete events (like neuron firing or 
exceeding a threshold of activation) occur at a state dependent rate 
frequently exhibit long-time persistent correlations, which are well 
reflected in power spectra of the representative signal recordings 
[13,15,30]. At the same time, many natural phenomena, including 
spontaneous human behavior, demonstrate non-homogeneous 
Poisson or even mixed, non-Poisson distributions of events, with 
intensity rates identified with sccde-invariant Levy statistics [31]. 

Levy-stable distributions are a class of self-similar infinitely- 
divisible probabihty laws [32,33] and by virtue of the Levy- 
Khinchin theorem are represented by distributions which can be 
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uniquely decomposed into two independent - Gaussian and 
Poisson - superposition parts. This structure is further responsible 
for the heavy-tailed amplitudal surges in the composition of Levy 



fluctuations and for long-ranged temporal dependencies charac- 
terized by the slow decay of their autocorrelation function. With 
this observation in mind, we have discussed the possible 
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emergence of a scale-invariant survival probability C{d) in a 
stochastic process of switching between two different activity states 
of the system (cf Appendix) pointing out that the two derived 
probabihty laws are related to each other through asymptotic 
properties of individual relaxation modes. 

Higher values of the exponent y for sleep-deprived subjects 
signals less heavy tails of waiting- time distributions in an immobile 
(resting) state than in an analogous distribution for the control 
group, and can be associated with restlessness/inquietude and 
increased variability (burstiness) of activity in recorded time series. 
Consequently, such alteration of locomotor behavior can be a 
representative sign of disorders related to sleep-deficiency and 
possibly, in line with results presented in former studies 
[8,16,17,34], a valuable diagnostic fingerprint discriminating 
between healthy and depressed/ disordered individuals. 

It seems likely that some aspects of both the waking period 
(length of prior wake, local use of neuronal networks) and the sleep 
period (sleep length and continuity) may explain performance 
deficits and the restoration of function. Accordingly, deeper 
theoretical und(;rstanding of different forms of sleep loss would be 
possible by comparing them together, in a single experiment [35]. 

Finally, it is important to recognize that the analysis of activity 
patterns of chronically sleep deprived subjects will not necessarily 
reveal their decreased capabilities. Compensatory efforts and 
recruitment of additional brain resources could intervene to make 
the artivity level appear luiafRxled. It could be that unraveling 
these disturbances requires more sophisticated analysis, or that the 
subtie neurocognitive effects of sleep loss can be observed only 
with the use of refined neuroimaging techniques, possibilities that 
deserve to be explored in future work. 

Appendix: Dynamics of Activation and 
Omnipresence of Power Laws - Theoretical 
Considerations 

We have studied the statistics of resting (inactivity) times and 
confirmed that the derived distributions are of power-law type 
C{d)cca^'' with the index y taking a fixed value over two orders of 
magnitude in time. In turn, activity period durations for control 
and sleep-deprived subjects obey a stretched exponential form for 
a wide range of recordings. We claim that both behaviors can be 
understood as different, asymptotic facets of an anomalous 
relaxation law governed by the generalized Mittag-Leffier 
probability distribution function describing rates of individual 
processes underlying switching from states of rest to activity. 

The pattern of time-duration (the overall pattern of distribution 
of fragments of time spent in a state of a given activity) is related to 
the rate function of a point process featuring events following 
spontaneous discharge (escape) from the inactivity (or alert) state. 
The spectrum of times between subsequent events can be 
described mathematically by a probability density function (PDF) 
pit) (which is equivalent to the experimentally derived frequency 
distribution of periods). For very small dt, tlu; probability that the 
event of escape from the state will happen in time range t to t + dt, 
given no such event occurred by times prior to t, can be expressed 
in terms of the conditional probability A(t)dt, where A{t) denotes 
the instantaneous "rate" at which the sequence of subsequent 
"ticks" of representati\'(; (■\'(;nts occurs in time. By introducing the 
probability that "the system remains in the state intact up to the 
time /", otherwise called the survival function C{t) 



C{t) = Pr{T>t)-- 



i: 



p{s)ds. 



(5) 



the conditional probability Pr{t< T<t + dt\T>t)= A(t)dt can be 
rewritten as 



A{t)dt-- 



l'+'"p{S)ds _ _ Pr{T>t+dt)-Pr{T>t) 



l^p(S)ds 



dPr{T>t) 
' Pr{T>t) 



Pr{T>t) 



dC{t) 



at) 



(6) 



This leads to a simple relation between the rate function A{t), 
the survival probability C{t), and the PDF p{t), which takes the 
form of 



A(0 = 



City 



(7) 



with the differential equation for the probability that no change of 
state occurs at least until time t: 



dCit) 
dt 



-Ait)Cit), 



(8) 



whose solution depends on the relevant choice of A(/): 

C(0 = Ce-Io'^«'''. (9) 



Note that cumulative probability that a state change happens by 
time t can be expressed as 



Fii) = Pr{T<t) 



(10) 



Accordingly, for te[0,oo) we will have C(0)=1 [i.e., C=l in 
Eq.(9)], lim,^oo C(0 = 0 and Cit)=\-Fit). From the above 
analysis it also follows that 



In particular, for a uniform rate A(?) = % = const > 0, one gets 



Cit) = e- 



Pit)-- 



[\-Fit)]=Xe- 



(12) 



which is typical for a homogeneous Poisson point process, where 
the probability A(/) = X is assumed constant throughout the whole 
time interval of interest. The mean duration time spent in a given 
state (say, before a discharge from the state takes place) is then 
given by Jg°° tpit)dt = 1 /A.. If, on the other hand, the time 
dependence of the rate function Ait) is assumed to scale dififerendy 
with small and large t as, e.g.: 



A(0 = 



akiXtf 



\+PiXtf 
which for 0<a< 1, p>Q leads to 



(13) 
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function. In this case the "effective" rate is derived from 



A(r)oc 



(Xtf 



the survival probability is given by 

City- ^ 



(14) 



corresponding to the frequency distribution of time durations 



(15) 



Pit)- 



a,l//i+l 



(16) 



The above formula Eq.(16) represents the probability density of 
the Burr distribution exhibiting two power-law behaviors 



d f " 



(22) 



In fact, it is easy to check that, if the rate PDF takes the Dirac- 
delta form 



(23) 



the relation Eq.(22) yi('lds a constant function A(t) = Xo and the 
corresponding PDF of duration times spent in a given state has an 
exponential form Eq.(12) with X = Xo. On the other hand, it can be 
shown that the rate Eq.(22) takes the form Eq.(14) if the intensity 
distribution function F/i{X) is of the generalized Mittag-Leffler 
form: 



iXtf-' : t^O it<<l/X) 
XO^cs ,„ , ■ (17) 

[(Xt)-"'!'-^ : t^oD it>>l/X) 



The survival function Eq.(15) possesses interesting limiting 
properties: As ^-»0, C(0-^e^*'''' , that is, it tends to the stretched 
exponential form (0<a<l) with the characteristic power-law 
behavior of the corresponding probability density for f< < l/X: 



p{t)az{Xt) 



,a-l 



(18) 



Note that as ;6->0 and a= 1, the common, exponential form of 
the survival function C{t) is recovered. In turn, for j8 # 0 and a = 1, 
the C{i) function tends to 



Fk{X)=Y,- 



(-1)T(/+-) 



=o/!r(-)r 



X 



a(/+l//i) 



, (24) 



which itself represents a mixture of the completely asymmetric 
Levy-stable SaiXifiyy^'") and the gamma gp(y)dy distributions 



where 



FdX)-- 



(25) 



(26) 



C(0- 



1 



[i+mxi)] 



(19) 



which yields a tail of the Pareto distribution in the expression for 
the probability density distribution p{t): 



p{i)oz(Xt) 



-l/H 



t>>l/X. 



(20) 



The nc)n-expc)n<;ntial form Eq.(15) and, consequently, the time 
dependent character of the rate A(?) may follow from a statistical 
probability mixing, e.g., from randomization of the rate X (cf 
Eq.(12)). 

To clarify this point, let us rewrite the survival function as a 
conditional probability 



PriT > t\A = X) = C(/| A = X) = e-^', 



(21) 



where A = X denotes random rate A taking value X. By taking an 
average with respect to PDF PAiX), the final C{t) takes the form 



C(t) = Pr(T>t} = ie-^'}= = f e-^'dFA(X), 

J 0 



Altogether, the above analysis classifies stretched exponential 
and Pareto distributions as probability laws stemming from the 
same origin, i.e.: either from the randomization of an individual 
rate X, or correspondingly, from the randomization of relaxation 
time r=l/X, the parameters which describe dynamics of 
stochastic activation events. The interpolation property of the 
Mittag-Leffler function (stretched-exponential character for short 
times with a transition to long-time inverse power-law behavior) is 
an appealing feature, and has been observed experimentally in 
various realms like protein conformation dynamics [36], dielectric 
relaxation in complex media [26,30,37] or financial market time 
series [38]. It well may be that the same function can serve in 
modeling temporal turnover effects in neurological recordings. 

Supporting Information 

File SI Details on data processing and fitting. The file 
includes analysis of goodness of fit of cumulative distributions for 
one and two activity thresholds; it contains 5 figures. 
(PDF) 

Acicnowledgments 

The authors express gratitude to Zdzislaw Burda for scientific discussions. 



i.e., it is given by the Laplace transform of the rate distribution 
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